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INTRODUCTION 


This  report  provides  a  simple  tool  for  the  orbital  mission  planner  to  use  in  designing 
optimal,  coplanar,  continuous- thrust  trajectories.  Six  charts  are  presented  that  relate 
vehicle  design  parameters  to  orbit  design  parameters  for  circle-to-circle  orbit  raising. 

Each  transfer  is  accomplished  with  constant  thrust  through  the  entire  flight  path  with 
no  restriction  on  its  magnitude.  The  direction  of  thrust  is  free  to  vary  within  the 
orbited  plane  with  fuel  expenditure  modeled  by  the  rocket  equation.  Intermediate 
eccentricity  is  allowed  to  grow  and  shrink  for  any  number  of  orbit  revolutions  as 
optimally  determined. 


There  has  been  a  substantial  amount  of  work  done  on  low-thrust  trajectories.  Of 
particular  relevance  to  this  report  are  the  works  covering  optimal,  continuous,  low- 
thrust  transfers  [Refs.  1  through  17];  included  in  these  are  numerous  papers  by 
Edelbaum  [Refs.  3,  7,  8,  and  11].  A  mission  planning  tool  for  optimal,  many-revolution, 
orbit  transfers  was  presented  by  Wiesel  and  Alfano  [Ref.  12]  where  an  analytically 
derived  solution  was  obtained  and  presented  in  graphical  form. 


The  equations  of  motion  used  in  this  optimal  control  problem  are  in  their  complete 
form  for  the  coplanar  circle-to-circle  case.  The  initied  acceleration  Aj  appears  in  these 
equations  and  is  allowed  to  vary  from  case  to  case  as  a  paraimeter.  Changing  the 
spacecraft’s  initial  acceleration  always  changes  the  minimum  time  to  accomplish  the 
orbit  raising.  To  present  the  solution  data  in  a  user-friendly,  universal  form  (Figs.  1 
through  7),  the  problem  is  rescaled  using  the  initial  orbit  radius  aind  gravitational 
parameter  of  the  central  body;  in  addition,  the  time-of-flight  t£  is  replaced  by  the  total 
accumulated  velocity  change  The  results  of  many  such  optimal  cases  are  presented 
in  graphical  form,  showing  the  relationship  between  spacecraft  acceleration,  propellant 
mass  fraction,  initial-to-final  orbit  radius  ratio,  and  the  minimum  accumulated  velocity  'or 
change.  The  graphs  themselves  are  composed  of  ptirametric  famihes  of  optimal  t 

solutions  assembled  in  comprehensive  charts  for  interpolation  by  the  orbitad  mission  ^ 

planner. 


The  charts  cover  the  tot£d  range  of  spacecraft  acceleration  values.  M«iny  users  of  space 
assets  are  interested  in  minimum-time  orbit  raising  and  repositioning,  so  the  high 
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acceleration  cases  will  be  of  increasing  future  interest  as  continuous- thrust  propulsion 
technology  improves. 


EQUATIONS  OF  MOTION 

The  differential  equations  that  define  spacecraft  motion  are  derived  under  the  following 
assumptions: 

—  The  force  of  thrust  is  constant  and  always  in  the  plane  of  motion. 

—  The  vehicle  has  a  fixed  propellant  mass  flow  rate. 

—  The  vehicle  acceleration  is  due  solely  to  the  force  of  thrust  and  a  spherically 
symmetric  inverse  square  central  gravitational  field. 

Using  the  defined  nomenclature  (page  13)  and  the  above  assumptions,  the  polar 
equations  of  motion  are 


f  =  u 


u  = 


V  = 


y!  _  £  X  ^i 
f  1  -f  m  t 

_UV  , 

r  1  -t-  mt 


(1) 

(2) 

(3) 

(4) 


where  the  dot  denotes  the  first  derivative  with  respect  to  time, 
rate  m  is  the  actual  mass  flow  rate  divided  by  the  initial  mass, 
when  dealing  with  propellant  loss. 


The  specific  mass  flow 
This  value  is  negative 


OPTIMAL  CONTROL  FORMULATION 

Given  an  initial  and  final  radius,  the  fuel  used  between  circular  coplanar  orbits  must  be 
minimized.  This  is  equivalent  to  minimizing  transfer  time  because  m  is  assumed 
constant.  From  Bryson  and  Ho  [Ref.  18],  the  Hamiltonian  can  be  written  as 

H  =  1  -I-  Af  f  -I-  Ay  U  -f  Ay  V  (5) 
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where  the  polar  angle  is  omitted  because  the  exact  final  position  is  not  specified  and  0 
does  not  appear  in  the  equations  of  motion.  The  behavior  of  the  Lagrange  multipliers  is 
given  by 


(6) 


Au  —  -  -^r  "b  V  (  r  )  (^) 

^v=  -Au(¥)+Av(^)  (8) 


The  partial  of  H  with  respect  to  the  control  parameter  ^  must  equal  zero 


m  _ 

d<j) 


^  cos  <i>  -  Ay  sin  ^ 


1  +  iht 


0 


(9) 


Knowing  the  acceleration  term  cannot  be  zero,  the  control  law  is  established  as 

tanc^  =  — ^  (10) 

-  Ay 

The  complete  transfer  is  characterized  by  the  initial  choice  of  Lagrange  multipliers.  To 
date,  the  only  closed-form  solution  that  exists  is  for  the  mainy-revolution  ceise  where 
intermediate  eccentricity  is  cissumed  to  be  zero  and  only  tangential  thrust  is  considered 
[Ref.  12].  Realizing  that  these  multipliers  can  be  scaled  without  affecting  the  control 
law  or  A  dynamics,  Aj.  is  set  to  -1  and  a  numerical  search  is  used  to  find  the  remaining 
two. 


ACCUMULATED  VELOCITY  CHANGE 

Accumulated  velocity  change,  V,  is  defined  here  as  the  total  velocity  impeirted  by  the 
force  of  thrust  during  the  elapsed  transfer  time. 


If 

=  f 

y  1  -I-  lilt 

0 


dt 


(11) 


This  formulation  of  the  well-known  rocket  equation  serves  two  purposes:  it  models  the 
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effect  of  fuel  depletion  and  also  defines  the  relationship  between  time  and  thrust- 
induced  velocity  change  through  the  equation 

A- 

“•^f  =  ^  In(l-l-mtf)  (12) 

Recasting  the  minimum-time  solution  in  terms  of  provides  the  means  for  a  very 
compact  graphical  representation  when  considering  the  initial  acceleration  Aj  and  the 
propellant  mass  fraction, 


nip  —  ”  rii  !/£ 


(13) 


Given  A-,  mp,  and  from  a  trzmsfer  chart,  the  specific  mass  flow  rate  aind  final  time 
axe  computed  as 


and 


In  ( 1  -  mp ) 


mp 

m 


(14) 

(15) 


For  the  limiting  case  where  m  and  mp  axe  zero,  the  totaJ  transfer  time  is  simply 


(16) 


CHART  GENERATION  AND  DISCUSSION 


Proper  scaling  can  eliminate  the  dependence  on  a  specific  central  attracting  body 
and  allow  a  global  mapping  of  solutions.  The  following  definitions  for  disteince  and  time 
units  axe  based  on  the  initial  radial  distance  r^  and  the  gravitational  parameter  of  the 
central  body  fi: 


1  DU*  =  rj 


1  TU* 


(17) 

(18) 
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Although  these  definitions  axe  dependent  on  the  physical  parameters  of  a  given  transfer, 
the  equations  of  motion  are  not.  Conveniently,  the  gravitational  parameter  is  always  1 
DU*3/tu*^  and  the  initial  values  of  the  (r,  u,  v)  array  for  any  circle-to-circle  coplcinar 
transfer  are  simply  (IDU*,  0DU*/TU*,  1DU*/TU*).  The  final  values  are  (RDU*, 
0DU*/TU*,  1/RDU*/TU*)  where  R  is  the  orbit  ratio  r£/rj-  For  this  study  only  orbit 
raising  is  considered  (r£  >  r- ). 

As  previously  mentioned,  the  initial  Lagrange  multiplier  values  completely  define  the 
transfer.  To  solve  this  two-point  boundary  value  problem  is  set  to  -1  and  a  shooting 
method  [Ref.  19]  is  used  to  find  A^,  Ay,  and  t£.  Minimum-time  solutions  for  a  variety  of 
orbit  ratios  and  mass  fractions  are  transformed  to  accumulated  velocity  change  md 
plotted  versus  initial  acceleration  to  produce  Figures  1  through  7  (pages  14  through  20). 

The  ripples  in  the  curves  of  Figures  1  through  4  reflect  a  transition  region  from 
gravitational  force  dominance  to  thrust  dominance;  they  can  be  used  to  define  the 
boundaries  of  low  thrust  (flat  curve),  intermediate  thrust  (rippled  curve),  and  high 
thrust  (constant  upward  sloping  curve).  Also,  the  ripples  and  subsequent  upward  slope 
testify  to  the  additional  cost  of  re-circularizing  the  final  orbit.  As  seen  in  Figure  5, 
fedlure  to  complete  a  revolution  during  transfer  increases  the  associated  cost.  The 
optimal  steering  law  causes  the  eccentricity  to  increase  for  the  first  half  of  each 
revolution  and  then  diminish  in  the  latter  half:  the  greater  the  thrust,  the  greater  the 
eccentricity  increase.  If  the  desired  radius  is  reached  before  completing  a  full 
revolution,  the  eccentricity  must  be  zeroed  out  to  meet  the  final  condition  of  a  circular 
orbit.  The  additional  cost  of  this  process  is  reflected  by  the  peaiks  in  the  curves.  When 
mahing  more  than  five  revolutions  of  the  central  body  this  effect  is  negligible  because 
the  induced  eccentricity  is  very  small  and  correctable  with  little  additional 
maneuvering.  If  the  final  state  is  reached  in  less  than  one  revolution,  the  intermediate 
eccentricity  and  associated  cost  grow  with  the  shortness  of  the  transfer  arc.  This  is 
reflected  in  the  positive  slope  of  the  curves.  The  cost  of  re-circulaxization  diminishes  as 
the  ratio  R  increases  because  the  gravitational  force  is  less,  making  vehicle  acceleration 
more  effective.  Also,  when  compared  to  Figure  1,  the  smaller  ripples  of  Figures  2 
through  4  reflect  the  increased  thrust  effectiveness  due  to  mass  loss. 

Figure  6  characterizes  all  transfers  where  Aj  <  10'“*  due  to  the  constancy  of  accumulated 
velocity  change  for  a  given  orbit  ratio;  this  is  in  agreement  with  the  results  of  Wiesel 


and  Alfano  [Ref.  12).  Figure  7  characterizes  all  transfers  where  Aj  >  4  due  to  the 
constant  slope  of  the  T£  lines  of  Figures  1  through  4;  the  intercept  value  Tq  is  used  to 
approximate  the  total  accumulated  velocity  change  through  the  empirical  equation 


(19) 


CHART  VERIFICATION  AND  USE 

As  a  matter  of  convention,  the  input  parameters  will  jJways  be  given  ais  the  array  (r^, 

Tf,  Aj,  mp).  It  is  assumed  the  reader  can  determine  total  accumulated  velocity  to 
two  significant  figures  from  the  transfer  charts;  six  significant  figures  will  be  carried  in 
all  computations  to  reduce  round-off  error.  Physical  constants  cam  be  found  in 
Reference  20.  The  data  will  be  presented  in  the  following  form: 

1)  Unsealed  Input  Array 

2)  Scaled  Input  Array  (Eqs.  17  and  18) 

3)  Accumulated  velocity  change  from  appropriate  chart  (Eq.  19,  if  needed) 

4)  Mass  flow  rate  (Eq.  14,  if  needed) 

5)  Minimum  time-of-flight  (Eq.  15  or  16). 

Test  Caae  ^1  The  first  test  caise  is  an  Earth-Mars  transfer  where  both  planetary  orbits 
are  assumed  to  be  circular  and  coplanar.  The  spacecraft  starts  in  a  heliocentric  orbit 
free  of  Earth's  gravitational  field  and  finishes  free  of  Mars’  field.  Earth  escape  and  Man's 
capture  are  not  considered. 

1) (  1.49598xl0"m,  2.27939x  10" m,  1.32712 x  10^° mVs^  8.33173 xlO’^m/s^,  0.0) 

2) {1.0DU^  1.52368DU*,  1.0DU*VTU*^  0.1405 DU7TU*^  0.0) 

3)  'V£  =  0.50  DU*/TU*  (From  Fig.  1) 

4) rii  =  0/TU* 

5) t£=3.55872TU*  =  1.78742x  lO^s. 

The  exact  numerical  answer  is  3.53186 TU*,  1.77393x lO^s,  or  205.316  days. 

Test  Case  ^2  An  Earth- Mars  transfer  solution  can  be  found  in  Bryson  amd  Ho  [Ref. 

18],  but  a  mass  flow  rate  of  -12.9  Ib/day  is  given  for  a  10“*  lb  vehicle  instead  of  the 
propellant  mass  fraction;  interpolation  will  be  needed  to  find  a  solution.  The  specific 
mass  flow  rate  is  -1.29x10  ^/day  or  - 1.49306 x  10’ ’^/s;  ais  a  starting  point  mp  is  set  at  0.25. 
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1)  ( 1. 49598 X 10" m,  2.27939x10" m,  1.32712 xl02°mVs^  8.33173 x  10'“ m/s^  0.25) 

2) (1.0DU^  1.52368  DU*,  1.0DU*VTU**,  0.1405  DU7TU*^  0.25) 

3) rj=0.54DU*/TU*  (From  Fig.  2) 

4) m=  -0.0748506  /TU*  =  -  1.49026x10' Vs 

5) t£=  3.33999TU*  =  1.67756x  10^  s. 

This  final  time  solution  is  adequate  because  rii  from  Equation  14  is  within  two 
significant  figures  of  the  original.  For  illustration  purposes,  mp  will  be  increased  to  0.5 
ajid  interpolation  will  be  used  to  refine  the  results  and  better  determine  the  transfer 
time. 

1) (  1.49598 X 10" m,  2.27939x  10" m,  1.32712 x  10^° m^s^,  8.33173 x  10'“ m/s^  0.5) 

2) (1.0DU*,  1.52368DU*,  l.ODU^VTU*^,  0.1405 DU*/TU7  0.5) 

3) r^=0.59DU*/TU*  (From  Fig.  3) 

4) rii=  -  0.165063 /TU*  =  -3.28637x10' Vs 

5)  t^=  3.02915 TU*  =  1.52144x10^  s. 

Linear  interpolation  of  rh  produces  a  final  ti  ae  of  3.33951  TU*,  1.67729x  lO^s,  or 
194.131  days.  This  is  within  two  significant  figures  of  the  Bryson  and  Ho  solution  of 
193  days  [Ref.  18]  and  the  exact  numerical  solution  of  192.748  days. 

Test  Case  This  case  involves  a  LEO-GEO  transfer  using  earth  canonical  units  and 
metric  units  for  a  propellant  mass  fraction  of  0.463,  requiring  interpolation  between  mp 
values  of  0.25  and  0.5. 

1) (1.O5DU0,  6.61DU0,  l.ODU®7TU0^  4.0x10'^  m/s^  0.25) 

2)  (1.0 DU*,  6.29524 DU*,  1.0DU*VTU*^  4.50079x10"'  DU*/TU*^  0.25) 

3) rj  =  0.60DU*/TU*  (From  Fig.  2) 

4) m=  -2.15799x1  O' VTU*  =  -2  48596xlO'Vs 

5) tf=  1.25849x10^  TU*  =  1.00565x  10®  s. 

Increasing  mp  to  0.5  yields 

1) (1.05DU®,  6.61DU®,  1.0DU®VTU®^  4.0x10'^  m/s^  0.5) 

2) (1.0DU*,  6.29524DU*,  1.0DU*VTU*^  4.50079x10'“  DU*/TU**,  0.5) 

3) r^  =  0.60DU*/TU*  (From  Fig.  3) 

4) m=  -  5.19952  X 10'“ /TU*  =  -  5.98972  x  10' Vs 

5)  7=  9.61628  X 10^  TU*  =  8.34763x10*  s. 

Linear  interpolation  to  match  the  given  mp  produces  a  final  time  of  1.00556 x  10^ TU*, 
8.72902x10*3,  or  10.103  days.  This  solution  is  close  to  the  10.0  day  result  of  Trussing 
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[Ref.  17]  and  the  exact  numerical  solution  of  10.0808  days. 

Test  Case  ^4  This  is  a  LEO-GEO  transfer  similar  to  ^3  with  the  acceleration  reduced 
by  a  factor  of  1000  and  a  propell^ult  mass  fraction  of  0.25. 

1)  (1.05 DU®,  6.61  DU®,  1.0DU®VTU®2,  4.0x10'®  m/s\  0.25) 

2)  (1.0  DU*,  6.29524 DU*,  1.0DU*VTU*2,  4.50079x10'^  DU*/TU*2,  0.25) 

3) r^=0.60DU*/TU*  (From  Fig.  6) 

4) m=  -  2.15799  X 10' VTU*  =  -2.48596x10'*% 

5) tj=  1.25849 X 10®  TU*  =  1.00565x10®  s. 

This  case  and  the  previous  one  involve  low  thrust  (flat  curve  for  Figures  1  through  4, 
constant  for  a  given  orbit  ratio  R  regcirdless  of  mp  or  A-).  Scaling  the  acceleration 
causes  cin  inverse  seeding  of  transfer  time.  As  expected,  the  final  time  is  1000  times 
greater  than  the  test  case  results  for  mp  =  0.25.  This  agrees  with  Reference  12. 

Test  Ca.se  ^5  This  is  a  high-thrust,  LEO-GEO  transfer;  it  is  similar  to  ^3  with  the 
acceleration  increased  by  a  factor  of  10®  and  a  propellant  mass  fraction  of  0.75. 

1) (1.05DU®,  6.61DU®,  1.0DU®VTU®2,  4.0xl02m/s2,  0.75) 

2)  (1.0 DU*,  6.29524 DU*,  1.0DU*VTU*^  45.0079 DU*/TU*^  0.75) 

3)  To  -6.4DU*/TU*  (From  Fig.  7),  rj=  42.9363  DU*/TU*  (Eq.  19) 

4) ih=  -  1.45309 /TU*  =  -  1.67292x10' Vs 

5) tj=  0.516144  TU*  =  448.05  s. 

The  exact  numerical  solution  is  445.582  s. 

Test  Case  ^6  The  last  test  case  is  an  intermediate-thrust,  Earth-Jupiter  transfer  where 
both  planetary  orbits  are  assumed  to  be  circular  and  copltinar.  Earth  escape  and 
Jupiter  capture  are  not  considered. 

1) (  1.49598xl0**m,  7.78299xl0"m,  1.32712xlO^®mVs^  1.77902 x  10'“ m/s^  0.75) 

2) (1.0DU*,  5.20260DU*,  1.0DU*7TU*^  0.03DU*/TU*^  0.75) 

3) rj=l.lDU*/TU*  (From  Fig.  4) 

4)  m  =  -0.037808  /TU*  =  -1.89897x  10®  /s 

5)  t£=  19.8371  TU*  =  9.96348xl0’'s  =  1153.18  days. 

The  exact  numerical  solution  is  1150.34  days. 
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CLOSING  REMARKS 


This  report  outlines  a  method  to  generate  minimum-fuel  circle-to-circle  coplanar 
transfer  trajectories  for  a  vehicle  with  constant  thrust.  Accumulated  velocity  change 
replaces  the  time-of- flight  while  also  accounting  for  propellant  mass  loss.  The  solutions 
are  globally  mapped  with  no  restrictions  on  initial  thrust  magnitude,  intermediate 
eccentricity,  or  number  of  revolutions  of  the  central  body.  Several  examples  axe 
presented  that  verify  the  transfer  charts  and  show  their  ease  of  use.  These  charts  are 
useful  tools  for  mission  planners  and  satellite  builders  to  assess  preliminary  fuel 
requirements  or  to  compare  different  propulsion  technologies. 
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NOMENCLATURE 

A-  =  initial  vehicle  acceleration 
DU*  =  scaled  distance  unit 
DU  0  =  earth  distance  unit 
H  =  Hamiltonian 

m  =  meters 

ih  =  specific  propellant  mass  flow  rate  of  vehicle 
J^p  ==  propellant  mass  fraction 

r  =  radial  distance  of  vehicle  from  attracting  center 
r£  =  final  radial  distance 
rj  =  initial  radial  distance 
R  =  ratio  of  final  to  initial  radial  distance 
s  =  seconds 

t  =  time 

t£  =  total  time  of  transfer  (time-of- flight) 

TU*  =  scaled  time  unit 

TU  0  =  earth  time  unit 

u  =  radial  component  of  velocity 

u  =  radial  component  of  velocity 

V  =  transverse  component  of  velocity 

T£  =  total  accumulated  velocity  chauige 

To  =  accumulated  velocity  change  intercept  value 

6  =  polar  angle 

\  =  Lagrange  multiplier 

Aj.  =  Lagrange  multiplier  associated  with  r 

Au  =  Lagrange  multiplier  associated  with  u 

Ay  =  Lagrange  multiplier  associated  with  v 

p  =  gravitational  constant  of  attracting  center 

4>  =  in-plane  thrust  direction  angle  (control  paxameter) 
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FIG.  1  Tf  contours  for  various  orbit  ratios  R  (mp  =  0) 
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FIG.  3  f,  coulours  for  various  orbit  ratios  R  (mp  =  0.5) 
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FIG.  4  fr  contours  for  various  orbit  ratios  R  (tnp  =  0.75) 
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FIG.  6  r,  contour  for  low  thrust  ( Aj  <  10 
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FIG.  7  Tq  contours  for  high  thrust  (Aj  >  4  ) 


